Polynomials of the best uniform approximation to 
sgn(x) on two intervals 



Alexandre Eremenko* and Peter Yuditskii^ 
August 24, 2010 

Abstract 

We describe polynomials of the best uniform approximation to 
sgn(x) on the union of two intervals [~A, —1] U [1, -B] in terms of spe- 
cial conformal mappings. This permits us to find the exact asymptotic 
behavior of the error of this approximation. 
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1 Introduction 

In [5] we obtained precise asymptotics of the error of the best polynomial 
approximation of sgn(a;) on two symmetric intervals [—A, —1] U [1, A]. Paper 

contains a somewhat simplified proof, together with generalizations. In 
this paper, we generalize the result to the case of two arbitrary intervals, 
the problem proposed to us by W. Hayman and H. Stahl, whom we thank. 

Related problems on the asymptotics of the error of the best uniform 
approximation by polynomials of degree at most n of the functions x"^^^ 
and l/(x — c), c ^ / on the union I of two intervals were completely solved 
by N. I. Akhiezer in [2]. 

Fuchs [3 [71 [8] studied general problems of uniform polynomial approxi- 
mation of piecewise analytic functions on finite systems of intervals. For the 
case of sgn(x) on two intervals I = [—A, —1] U [1, B], the result in |6] is 
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Here 



Ln = inf sup |sgn(rc) — 



where Vn is the set of polynomials of degree at most n; positive constants 
Ci and C2 depend on A, B, and r/ is the critical value of the Green function 
G of the region C\/ with pole at infinity. The arguments in [6] do not give 
optimal values of Ci, C2. 

When A = B,we have = {A — 1)/{A + 1), and the result obtained 
in [5] is 

W = Wi - ^^'^ 2m + l^^y . (1.1) 

In this paper we will obtain a result of the same precision for arbitrary A 
and B. In the case A ^ B, the ratio ^/ne'^'^Ln oscillates. Similar oscillating 
asymptotic behavior was found by Akhiezer for the polynomials of least 
deviation from 0, that is, for the error of the best uniform approximation of 
j-n+i |-|y polynomials of degree at most n on two intervals. 

To state our main asymptotic result we introduce certain characteristics 
of the region C\/. Let 

/X C — X 

— dx, — 1 < X < 1, 

-1 ^{l-x^){x + A){B-x) 

be the Green function of C\/ with pole at infinity, see, for example, [3], 
where C G (—1, 1) is the unique critical point, 

_ - + A){B - x))-^l'^xdx 

~ - x2)(x + A){B - x))-V2dx ■ 

We introduce positive constants r] = G{C, 00) and 

m = -Ig"{C) ^ 



2 2^{l -C^){G + A){B-Gy 

The Green function G{z, C) satisfies 

G{z,C) = -\n\z-G\+r]2 + 0{z-C), z^G, 
and this relation defines the Robin constant 772- 
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Let uj{x) = uj{x, [—A, —1], C\/) be the harmonic measure of the interval 
[—A, —1]. An exphcit formula for to is 



uj{z) = Im 



f\{{x^ - + A){B - x))-i/2dx' 
In our notation related to theta functions we follow Akhiezer's book [3]. 



(1.2) 



Theorem 1.1. The error Ln of the best polynomial approximation o/sgn(x) 
on I = [— A, — 1] U [1, B] satisfies 



Ln = (c + o(l))n" 



-l/2g-nr, 



^o(^({na;(oo)+a;(C)}-a;(C))|T) 
i9o^({na;(oo)+a;(C)} + a;(C))|r 



(1.3) 



where 



c = 2(^r/i)-i/2e-'?2^ 



jfi(t'-l)iB-t)iA + t))~y'dt 
''j\i{l-t^)iB-t)iA + t))-y^dt' 



(1.4) 



and 



'do{t\T) = l- 2h cos 27rt + 2/1"^ cos Airt - 2h^ cos 67rt + 



h 



is the theta- function. In (1.3) we used the notation {x} for the fractional 
part of X. 

Our method is somewhat different from the methods of previous au- 
thors. It is based on an exact representation of the extremal polynomial as 
a composition of conformal maps of explicitly described regions. This can 
be considered as a development of the arguments in [U [5] . Our representa- 
tion of extremal polynomials permits to find their asymptotic behavior in 
various regimes and their zero distribution. Actually, the main asymptotic 
result of this paper is Theorem 7.1, which has somewhat technical state- 
ment, and Theorem 1.1 is a simple corollary. For example, according to [S], 
the numbers ni and n2 of zeros of the extremal polynomial Pn on [—A, —1] 
and [1,-B], respectively, satisfy 

lim ni/n = uj{oo) and lim n2/n = 1 — a;(oo), 

while our Theorem 7.1 implies a stronger conclusion: ui = na;(oo) + 0(1). 

However, in this paper we focus on the error term of the polynomial 
approximation, and do not explore other corollaries from Theorem 7.1. 
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Figure 1: Graph of the function S{z,a) on the real axis 



A representation of extremal polynomials is described in sections 2, 3, 
where we use an entire function introduced in [5]. In Section 4 we find an 
integral representation of the principal conformal map involved, and then 



study its asymptotics in sections 5-7. We derive (1.1) as a special case of 
Theorem 1.1 in Section 8. Finally, in Section 9, we sketch without proof the 
limit case B = oo. In this case, instead of approximation by polynomials 
one has to consider approximation by entire functions of order 1/2, normal 
type. 



2 Preliminaries 

We begin by recalling the construction of the entire function S{z, a) of ex- 
ponential type one which gives the best uniform approximation of sgn(x) on 
the set (— oo, —a] U [a, oo), where a > 0. There is a unique such function for 
every a > 0; it is odd and satisfies 

S{a, a) = 1 - L(a), 5(4, a) = 1 - (-l)^L(a), (2.1) 

where L{a) is the approximation error, and ci < C2 < . . . the sequence of 
positive critical points. The graph of this function is shown in Fig. 1. We 
define the positive number b = b{a) by cosh 6 = 1/L(a). It is easy to see that 
6 is a continuous increasing function of a, and the correspondence a i— )• 6 is 
a homeomorphism of the positive ray onto itself. 
For every 6 > 0, we consider the region 

= {x + iy : X > 0,y > 0,x > arccos(cosh6/ coshy) for y > b}. 



4 



71/2 n 2n 3n 



Figure 2: Domain such that S{z) = 1 — L{a) cos'i/'(z), -0 '■ 



n. 



This region is shown in Fig. 2; it consists of the points in the first quadrant 
to the right of the curve 
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{arccos(cosh 5/ cosht) + it : b < t < oo}. 



Let ip be the conformal map of the first quadrant C+_|- onto 0, normahzed 

by 

tpiz) = z + z oo and ■0(0) = b. (2.2) 

Put a = V^-HO). 

In [5j we proved that 

S {z, a) = 1- L{a) cos i){z), z G C++ := {z : Re z > 0, Im z > 0}. (2.3) 



As the right hand side of (2.3) takes real values on the positive ray and 



imaginary values on the positive imaginary ray, S extends to the whole 
plane as an odd entire function. 
The following asymptotics hold 



lim \/ae°'L{a) 



(2.4) 



Notice that all critical points {±Cfc} of S{z^ a) are real, and ip{ck) = irk. 
It is convenient to modify a little this conformal mapping. We write 
S{z,a) := S {Vz^~+~a? , a) , where z belongs to the upper half-plane C+ with 
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Figure 3: Preimage S (M, a) in the upper half-plane 



the slit {it : < t < a}, see Fig. |3j Function S is not entire, it is only 
defined in the upper half-plane. 

Again we have the conformal mapping ^ : C++ — )■ Q but in contrast 



with (2.2) 



'iP{z) = z + ..., z^oo, V(0) = 0, (2.5) 

and therefore ■ip{ia) = ib. The full preimage of the real axis under S in the 
upper half-plane consists of the curves 



4 = ^-i({7rA; + it:t> 0}), k = ±l,±2,.... 



(2.6) 



shown in Fig. 3. These curves have vertical asymptotes {Re z = irk — 7r/2}. 

Now let Pn{z) be the best approximation of sgn(j;) by polynomials of 
degree at most n on two intervals / = /_ U /+, /± C M±. Using a linear 
transformation we may always assume that / = [—A, —1] U 

Our goal is to obtain a representation for the extremal polynomial in the 
form of the composition 

P„(z) = 5(e„(z),a„) (2.7) 

where 0„ is the conformal mapping of the upper half-plane on a suitable 

^In what follows, the letters O and 8 are used to denote conformal maps which have 
no relation to theta-functions i?. 
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Figure 4: Domain 112,3 



"curved" comb-like region, and a„ is an appropriate value of the parameter a. 



First we give typical examples of the representation (2.7) and then show 
that these examples exhaust all possibilities. 

First example. For n = 4, consider the following region 112,3, see Fig. |4[ 
Its boundary consists of the vertical segment [0,ia], the horizontal segment 
[— C2,C3] and the curves 8-2,63. 

Let Q{z) = 64(2;) be the conformal mapping 

e : C+ = {z G C, Im z > 0} ^ n2,3, e(±l) = 0, G(oo) = 00. (2.8) 

The function P{z) = S{@{z), a) can be extended to the lower half-plane due 
to the symmetry principle. Therefore it is an entire function, which is, in 
fact, a polynomial of degree 4 due to its asymptotics at infinity. The graph 
of this polynomial on the real axis is of the form given in Fig. [5| where 
(j) = By the Chebyshev theorem (for the two interval version of this 

theorem see [1], [2] [3]), P{z) is the extremal polynomial on the set / with 
A = -0(-C2), and B = 0(03). 

Second example. Let us point out that the above polynomial P{z) has 7 
points of alternance, instead of 6, which are required by the Chebyshev 
theorem for a polynomial of degree 4. Therefore the same polynomial is 
extremal on the sets of two kinds 

I =[- A, -1]U[1, (Pics)], (/>(-C2) < < 0(-ci) (2.9) 



or 



-C2),-1]U[1,5], 0(C2) < 5 < 0(C3). (2.10) 
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Figure 5: Extremal polynomial corresponding to the region 112,3 



Third example. From the position / = [(/){— 02), —1] U [1, 4'{c2)] we can start 
a deformation of the set / and of the extremal polynomial. Namely consider 
the region n^3(/i), see Fig. 6 Here we added to the boundary a segment of 
the curve 62 that starts at the critical point C2 and has length h G (0, 00). 
In this case 

e : C+ = {z G C, Im 2 > 0} ^ ^tsW^ ©(±1) = 0' ©(oo) = (2-11) 

and P{z) := S{@{z),a) is of the form given in Fig. [7| For this new family 
of regions, which we parametrized by positive h, the polynomial is extremal 
on the set I = [(/)(-C2), -1] U [l,(/)(c2 - 0)]. 

In the next section we show that these examples exhaust all possibilities 
for the extremal polynomials. 

3 Parametrization 

We begin with a general description of extremal polynomials. Fix a > 0. 
This defines the number L = L{a) and the function S{-,a). Let ki and k2 
be two positive integers. Consider the region Ilki,k2 the upper half-plane 
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Figure 6: Domain n^3(/i) 
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Figure 7: ...and the corresponding extremal polynomial 
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bounded by the curves and (5^2- Then for /i > and k2 > 2, we define 
the region (h) by making in I^kiM ^ ^^^^ along 6^2-1 starting from 0^2-1 

and such that the length of its image under ip is h. So n^_^ ^3(0) = ^ki,k2- 
Similarly we define 11^^ fcj^^^ ^1 > 2 by making a slit along 

Let be the conformal map of the upper half-plane onto Il^_^j^^(h), 
normalized by G(±l) = 0, 6(00) = 00. Consider the function 

Piz) = S{e{z),a). (3.1) 

By construction, it is real on the real line, so the symmetry principle implies 
that P extends to an entire function. By looking at the asymptotic behavior 
as 2; — )• cxD we conclude that P is a polynomial of degree fci + /c2 — 1- All 
critical points of this polynomial are real. If /i = 0, then all critical values 
are — litL and libL on the negative and positive rays respectively. If /i > 0, 
the extreme left critical value is changed to —1 it Lcosh/i, or the extreme 
right critical value is changed to 1 it L cosh h. 

We have seen in the previous section that each of these polynomials P 
is the extremal polynomial for some A and B. Now we prove that for every 
A and B one of these polynomials is extremal. 



Proposition 3.1. All extremal polynomials are of the form ( 3.1) with Q as 
defined above and some ki,k2,a and h. 

We give an elementary proof of this proposition, which is based on count- 
ing critical points and alternance points. This proof does not extend to the 
case of entire functions, so in Section 9 we will give another proof which is 
less elementary but avoids counting. 

Proof. In the proof we will use the following fact which is well-known and 
easy to prove. 

Let Pi and P2 be two real polynomials with all critical points real and 
simple, and suppose that their critical points are listed in increasing order 
as ci < C2 < ■ ■ . < Cn~i and c[ < C2 < . . ■ < c'l-i- V Pi{cj) = P2{c'j) for 
1 < J < — 1; then Pi{z) = P2{cz + b) for some c > and real h. 

For a discussion and generalizations of this fact to entire functions, see 

Let A> l^B > \, and a positive integer n be given. (We will deal with 
the degenerate case A = lox B = \ later) . Let P be the extremal polynomial 
of degree n which exists and is unique by Chebyshev's theorem. According 
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to Chebyshev's "alternance theorem" , this polynomial P is characterized by 
the following properties: let Q{x) = P{x) — sgn(x), then 



|Q(x)|<L, xG [-A-1]U[1,5], (3.2) 



and there exist 



m>n + 2 (3.3) 
points xi < X2 < . • • < Xm in [— ^) —1] U [1, -B] such that 

|Q(xj)| = I < j < rn, and Q{xj)Q{xj+i) < 0, I < j < m — 1. 

(3.4) 

These points xj are called the alternance points. Evidently, all alternance 
points in {—A, —1) U (1, B) are critical, that is, P'{x) = at all such points. 
Let K be the number of critical alternance points and the number of 
non-critical alternance points. We have the evident inequalities 

K <n-l and < 4. 



Combined with (3.3) this gives 

n + 2<m = K + N<n + 5. 

So we have three possibilities: 

a) m = n-|-3, A^ = 4, K = n — 1. The last two equalities imply that all 
critical points of P are real and simple, and each of them is an alternance 
point which belongs to (—^,—1) U (1,-B). All 4 points —A, —1,1, B are 
non-critical alternance points. So the graph has the shape shown in Fig. 5. 

b) m = n -|- 2, = 3, K = n — 1. Again all critical points are real, 
simple, belong to [—A, — 1] U and each critical point is an alternance 
point. All endpoints —A, —1, 1, B except possibly one are alternance points. 
Let us show that —1 and 1 are alternance points. 

Proving this by contradiction, suppose, for example that —1 is not an 
alternance point. Then 1 cannot be a critical point because = 3. Thus 
P'{x) / on [-1, 1] and P(l) >1-L> -1 + L> P{-1), we conclude that 
P is strictly increasing on an interval (—1 — e, 1 -|- e) for some e > 0. This 
implies that 1 is also not an alternance point, a contradiction. 



Thus the polynomial P is of the type described in (2.9), (2.10). 

c) m = Ti + 2, K = n — 2, N = 4. In this case we have exactly one simple 
critical point z which is not an alternance point. Evidently this exceptional 
critical point is real. We claim that it belongs to M\[— A, B]. 

First of all, z ^ {—A, —1,1, B} because A^ = 4 so all these 4 points 
are non-critical. Second, z cannot be in the interior of one of the intervals 
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{—A, —1) or {1,B). Indeed, if it is in the interior of one of these intervals, 
consider the adjacent alternance points xj and Xj+i on the same interval 
such that Xj < z < Xj+i. Such Xj and Xj+i exist because all endpoints of 
each interval are alternance points, and z is not an alternance point. As z 
is the unique critical point on {xj,Xj+i), we obtain a contradiction with the 



alternance condition (3.4). Finally we prove that z ^ (—1,1), Proving this 
by contradiction, suppose that z G (— 1, 1). As —1 is an alternance point we 
have P(— 1) = — libL. Suppose first that 

P(-l) = -l-L. (3.5) 

Then P'(— 1) < because —1 is not critical (A^ = 4 in the case we consider 
now), and ( |3.2| ) implies that 1) < 0. As P' changes sign exactly once on 
(—1,1) and the point 1 is also non-critical, we conclude that -P'(l) > 0. As 



P(l) = 1±L, (3.2) implies P(l) = 1 — L. This equality and (3.5) contradict 



the alternance condition (3.4). The case P{—1) = — 1 + L is considered 
similarly. 

Let c be the critical point which is outside [vl, -B]. It is easy to see that 
|P(c) + 1| > L if c < and \P{c) - 1| > L if c> 0. This is because = 4 
and thus A and B are non-critical alternance points. 

So in the case c) we have the graph of the type shown in Fig 7. 

To summarize, we proved that in all cases the critical points are real 
and simple, all critical values, with at most one exception are — 1 it L on 
the negative ray and 1 it L on the positive ray, and the exceptional critical 
value, if it exists, corresponds to an extreme (left or right) critical point. If 
the exceptional critical point c is positive then |P(c) — 1| > L and if c is 
negative then |P(c) -|- 1| > L. 

Polynomials S{Q, a) constructed above permit to match any such critical 
value pattern, so we conclude that P{z) = S{Q{cz + b),a) with c > and 
6 E M. Finally, the points —1, 1 are always non-critical alternance points, 
and this implies that c = 1 and 6 = 0. 

It remains to consider the degenerate case. Suppose, for example that 
B = 1. Then only 3 alternance points can be non-critical, so we are in 
the case b). The extremal polynomial in this degenerate case can be easily 
written explicitly: 

where Tn{x) = cosrearccos a;, and 

Ln = 3— TTTlt n\ ~ 4exp ( — nch~^ ( 1 + 



r„(i + 4/(^-i)) + i "V V ^-1 
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is the approximation error. 



□ 



Remarks. It is easy to see that our polynomials depend continuously on 
h. When /i oo, we have 11+ ^^^(/i) Ilk^^ki-i and 11" n^^.i^fc^ 
in the sense of Caratheodory, and the corresponding polynomials converge 
uniformly on compact subsets of the plane. 

Let us show that A{h) and B(h) depend monotonically on h. Let h > hi, 
e{z) = e{z, h) and Gi(z) = G(z, hi). Then the function w{z) = o e{z) 
maps the upper half-plane into itself, and has the properties: tt)(ibl) = 
±1, w{co) = oo. Thus it has a representation 

w{z) = p^{z-l) + l + j _ ^^ dp{x) 

where dp is positive and supported on a compact set / such that x > B for 
all X £ I, and poo > 0. Here we used w{l) = 1 and 10(00) = 00. Now we use 
the condition w{—l) = —1. We obtain 

Poo = 1 - / + (3-6) 

Since w{B) = Bi we have 

Bi-B = {p^- 1){B - 1) + - 1) y 



Using (3.6) we get 



B -\- \ 

Bi-B = {B-1) I — — -dp{x) > 0. 

^ ' ' {x- l){x - B){x + 1) ' 



Similarly 



-Ai + A=iA + l)l > 0. 



4 Integral representations 

Asymptotic relations for the extremal polynomials are based on an integral 
representation of the conformal map 0. 

Consider the conformal map of C\I onto on the annulus in Fig. ^ Here 
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Figure 8: The confer mally equivalent annulus. 

we assume that the upper lialf-plane is mapped on the upper part of the 
annulus with the following boundary correspondence 

B ^ -1, -A ^ -A, -1 A, 1^1. (4.1) 

By G{z, zq) we denote the (real) Green function of the region C \ /, 
/= [— A, — 1] U [1, -B], with pole at Zq. In particular G(z) := G{z,oo). Recall 
that in the upper half-plane G{z) can be represented as the imaginary part 
of the conformal mapping ^(z) of the upper half-plane onto the region 11, 
Fig[9| 

G{z) = Im $(z), Im z > 0, $(±1) = 0, $(oo) = oo. (4.2) 

The map <&(z) defines certain important characteristics of the region: 
the critical value 

rj = G{C), C G [-1, 1] such that VG(C) = 0, (4.3) 

and the harmonic measure uj{z) of the interval [—A, —1]. We have 

^(—A) = — vra, where a = uj{oo). (4-4) 

Now we associate to ^{z) the function 

g{C) = -i^ziO) (4.5) 
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Figure 9: The image 11 of the map 



where C belongs to the upper half of the annulus, see Fig. [8] This function 
can be extended to the upper half-plane by the symmetry principle. We 
have G{z(()) = Re g{C)- We call g the complex Green function. 

Let —fi G [—1, —A] corresponds to the infinity in the z-plane, —/j. = C(oo) 
(see Fig. [s]). We define the jump function 

\0, «£(-(., -A) 

which we extend by the symmetry = j(A^^) = on the whole 

negative ray. Since 

r 0, e>o, 

[ 7r(a- 1), e G (-1,-/^), 

we obtain the following integral representation for g{(^) in the upper half- 
plane 
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Remark. In the representation (4.7) the normahzation condition ^{1) 
was used. The second normahzation condition <&(—!) = gives 



a = J 



._J-oo{g^-g^}j-m 



In what follows we will use the bar over a function to denote similar averages. 



Naturally we can simplify (4.7), but the point is that we can write a 
similar representation for the conformal mapping Qn{z). Recall that for a 



given n, there exists a unique region n(n) = n^^^.^^ see Fig. 

such that the conformal mapping G„ : C+ — )• n(n) represents the extremal 
polynomial (2.7). We define the function 

e„(C) = -ie„(z(c)), ^„(a2c) = 0„(C). 

We write the imaginary part of ^^(C)) < a sum 

vr 



(4.8) 



Im 9n{0 = T^hin) --+ TinjiO + Xn{0, 



(4.9) 



so that XniO is a continuous function, which is normalized by the condition 
Xn{-lj) = 0. Then 



.(C) 



1 



1 



{'Kki{n) 



vr 



Theorem 4.1. In the above notations 



T^nm+xnimi- (4.10) 



AO-ng{0 



1 



(Xn(0 - Xn)di, 



(4.11) 



where 



Xr. 



p = vrna — 7rKi(n) H — . 



(4.12) 



Proof. We subtract ng{C,) in the form ( |4.7| ) from (4.11). Then, we use the 
second normalization condition g[\) = 9n{\) =0. □ 

This representation will imply Fuchs' asymptotics as soon as we show 
that Xn(0 is uniformly bounded. 
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5 Fuchs' asymptotics 

Let us begin with a simple remark. 

Lemma 5.1. Let w{z) be a conformal mapping of the upper half-plane onto 
a sub-region of the upper half-plane which contains the half-plane Im w > tq. 
Assume the normalization w{z) ^ z, z ^ oo. Then 

Im w{z) — Im z £ (0, tq). (5-1) 

Proof. Consider the integral representation of Im w{z) 

1 r"^ 

Im w{z) =lm z-^ - P{z,t)v{t)dt, (5.2) 

where P{z, t) is the Poisson kernel. Since < v{t) < tq we obtain the desired 
inequality. □ 

Proposition 5.2. There are constants Ci and C2 (depending of the given 
system of intervals I ) such that 

Ci < On — nrj < C2. (5-3) 

Proof. Recall that the curves in Fig. [3] were defined as preimages of vertical 
lines in the region Q in Fig. [2] under a conformal mapping which maps the 
right half-plane into a sub-region of the right half-plane. Thus we can apply 
Lemma 5.1, to obtain that |Xn(OI ^ ^vr. Therefore |xn(0 — Xn\ is also less 
than 27r. 

Now, On is the maximum of ^n(0 o^i the interval (A, 1) and r] is the 
maximum of g{^) on the same interval. Due to the integral representation 
(4.11) the difference between these two functions is uniformly bounded in 
this interval. Thus (5.3) is proved. □ 

Corollary 5.3. The following limit relation holds 



lim^=5(C), (5.4) 
n— >-oo n 



in particular 



lim ^ = a, lim — ^^-^ — — = lim — = rj. (5-5) 

n— >oo n n— >oo n n— >oo n 



Proof. We divide (4.11) by n and pass to the limit. □ 



17 



Figure 10: The rescaled region n(n)/n for a large n. 



Corollary 5.3 has the following geometric interpretation. Making the 
rescaling n(n) — t- Il{n)/n we obtain the limit conformal mapping onto the 
region shown in Fig. |9j Let us look more carefully at the limit procedure, 
see Fig [To} the distance between the additional cut and one of the infinite 
cuts (left or right one) approaches zero, however the position dn of the end 
point of the additional cut influences the asymptotic behavior along various 
subsequences {ni}. We define the subsequences by the condition: there 
exists a limit d = d({n/}) = lim/_^oo dnr Taking into account this point d, 
in the next section we describe the asymptotic behavior of L„ in a more 
precise way. 



6 The limit density xiO 

We fixed a subsequence {ni} such that the limit d = lim/^oo dm exists. Our 
main goal in this section is to show that the limit density 



lim Xn, (0 

t— >oo 



(6.1) 



exists and to find this limit. 
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We start with the following general lemma. Let /, be a bounded increas- 
ing difFerentiable function defined for x > 0, and suppose that f{x) = for 
< X < 6 with some 6 > 0. We consider the region 

nf = {z = x + iy:x>0,y> f{x)}, (6.2) 

(it looks like the region $7 in Fig. [2] reflected in the line x = y). Let w be the 
conformal map from the first quadrant C++ onto / with the normalization 

w[z) ^ z ^ CO, 'w{0) = 0. (6.3) 

Let a be the point such that w{a) = b. 

Lemma 6.1. Let w{x) = u{x) + iv{x), x > a. Then f{x) < v{x). 

Proof. We extend w by the symmetry principle to the map of the upper 
half-plane into itself (the extended map is still denoted by w), and use the 
integral representation 

For X > a we have 

, , 1 2x v(t) — v(x) , v(x) , x + a . , . 
w{x)=x + - ^—^dt + -^\n +iv{x). (6.5) 

IT J a t + X t — X TT X — a 

Therefore 

1 f°° 2x v(t) - v(x) , v(x) , x + a , , 

u(x)=x + - — —dt + ^^ In > X. (6.6 

■K J t + X t — X TT X — a 

Since f{x) is increasing we obtain 

fix) < f{u{x)) = v{x). (6.7) 

□ 



We apply Lemma 6.1 to obtain the limit density for the conformal map of 
the first quadrant onto the region Vt in Fig. 2. Namely, as before we consider 
the conformal map w{z) = —iip{—iz), where ■0 is defined in (2.5) and ex- 
tended by symmetry to the right half-plane, and the integral representation 
(6.4) for it. Let us notice that in our case we have the exact formula 

f{x) = arccos ^ , x >b. (6.8) 
coshx 
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Between the values a and b there is a one-to-one correspondence, moreover 
b ~ a + 1/2 In a. Thus we have the density v{x) = v{x,a) in (6.4) as 
a function of the parameter a and we are interested in the hmit density 
v{x) := hma^oo v{ax, a). 



Lemma 6.2. The following limit exists 

lim v{ax, a) 



V(X} 



0, X < 1 
f, x>l. 



(6.9) 



Proof. It is evident, that v{x) = for x G (0, 1). For x > 1 we use Lemma 6.1 
and the asymptotic relation between a and b: 



v(ax, a) > arccos 



cosh b 



arccos ■ 



cosh ax g(x-i)a 
On the other hand v(ax, a) < vr/2, thus the lemma is proved. 



(6.10) 

□ 



Now we are in position to evaluate the limit density (6.1). 



Theorem 6.3. Let {n^} be a subsequence such that limdn^ = d. Without 
loss of generality, we assume that He d = 7r(l— a) (alternatively Yie d = ira). 
The relation 

d = <^{D)=ig{-K), 



uniquely defines D G [-B, oo] and —k G Then 



where y = Re g{£,). 



vra 



2 7|t|<,, {t - vY + (vra)' 
1 f 7r(a-l) 



-,dt, 



'^J\t\<^ (t-y)2 + (7r(a-l))2 
1 f 7r(a-l) 



[2 7|,|<^ {t-yf + {7:{a-l)f 



dt, 



dt + TT, 



(6.11) 

-/u < ^ < -A, 

-K < ^ < — /i, 
-1 < ^ < -K, 

(6.12) 



Proof. First we assume that — ^ < C < L^t zi = (^). For a sufficiently 
large /, by (6.4), we have 

If Im zi 



Im wi = Im zi + 



vr J\t\>a„. (t - Re z;)2 + Im z 



^v{t,ani)dt. 



(6.13) 
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Figure 11: The preimage of the level line Im w = t 



(see Fig. 11) we obtain 
1 



Substituting Im zi = TTki{ni) — f + Xni(C) aiid I™ = '^ki{ni) to (6.13) 

(7rfci(n/) -Tr/2 + Xni{0)'"{t,ani) 



vr/2-Xn,(e) 



1 /■ {irkiim) -1T/2 + XniiO)0'niVianit,ani) 



^ ^|>i - yii(.W + (7rA;i(ni) - 7r/2 + XmiO) 



dt 
dt 



By the leading term asymptotics, Corollary |5.3[ we have 

ki{ni) 



lim 

;-5-oo ni 



a, lim — = r/, lim = y\c,)- 



Passing to the limit in (6.14) we get 
vr/2 - x(6 ^ 



rjTTa 



;V{t)dt, 



(6.14) 



(6.15) 



(6.16) 



By Lemma 6.2, after trivial simplifications we obtain the first equation in 



(6.12). 



In the second case —k < ^ < — /i, for sufficiently large I, the point zi 
corresponds to a point wi on the line Im wi = TTk2{ni). Thus we can repeat 
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the previous arguments with a replaced by 1 — a (let us mention that Xni (0 
is negative here). 

In the last case Im wi = 7r{k2{ni) — 1), and this leads to the shift of the 
limit value by tt. □ 



7 Simplifying the result 

In this section we prove the following theorem, which is our main result, and 
which implies Theorem 1.1. 

Theorem 7.1. Let v he the 'point in the interval (A, 1), such that g{v) = rj. 
Fix a subsequence {n/} such that lim/_^oo dm = d = ig{—K). Let g{C, v) and 
g{C, —k) he the corresponding complex Green functions. Then 



1 - ff(C) 



hm {Om (C) - nwiO] = - In + <?((, y) " 5(C, -'^). (7.1) 



Proof. First of all we split x{C) into the sum of a continuous function Xc{C) 
and the jump 

f 1, f G (-1, -k) , , 

[0, ^G(-K, -A) 

As usual the jump function is extended to the negative ray by the reflections 

ji(i/6=ii(0,Ji(A'e)=ii(0- 

Note that the jump function is related to the Green function G{z,D), 



compare (4.6) and (4.7). Since 



0, e > 0, 

Im ^(C, -k) = <( TTUj{D), ^G(-k, -A), 

7Tiu{D)-l), ^e(-l,-K), 

we obtain the following integral representation for g{(^, —k) in the upper 
half- plane 

^(C,-'^) = - ^} Cji-JiiO)dt (7.3) 

where we use the notation j'l introduced in the Remark in Section 4. Notice 
that G{z{C),D) = Re ^(C, -k), and 

Ji = uj{D). (7.4) 
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gO) 





gH) 



Figure 12: The image of the Green function 



Due to the chosen normahzation xi~f^) = 0) '^^ have 

The main point is to evaluate the Cauchy transform of the continuous 
part Xc(0- 

Lemma 7.2. Let v G (A, 1) 6e such that g{v) = rj, that is, v corresponds 
to the critical point C G (—1,1). Let giC^v) he the corresponding complex 
Green function. Then 

^ ^ ^ ^ ^ {xM)-Xc)di=\\n'l^^^+g{(:,v). (7.6) 



^7-00 U-C e-lj'" '"' " 2 l^ + g{C) 



Proof. Using (6.12), for z = y{^) + ia = g{^), ^ G (— —A), we have 



So let us consider the function 



/(C):=^ln^^. (7. 

2 -ri-9{C) 
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Figure 13: The image of {rj — g{C))/{—V — oiO) 



The image of the function g{C) is shown in Fig. 12 the image of the fraction 
hnear transformation is shown in Fig. 13 Let us point out that for C in 
the upper half of the ring, Fig. |8| we obtain the values of in the right 
71 — giC) 

half-plane and for 7^ in the unit disk. Thus, 



piO ■■= Im fiO 



7r/2, A<C<z^ 
-7r/2, u <^<l 



(7.9) 



here i' is such that g{L') = rj. 

We use the integral representation of f{() + i7r/2 



1 , v-g{0 ft.. , TT. 1 



1 



1 



vr 



P(0 + 2 j 



(7.10) 

The complex Green function related to the critical point C G (—1,1) we 
still normalized by the condition g[l,v) = 0. Therefore 

-vr, ^ G (A, I/), 

luig{i,v)={ -vr(l-u;(C)), ^ < 0, 
0, 



According to (7.9) we can represent g{C,v) as 



1 



e-c e-i 



1 



1 



e-c e-i 



vr 



P(0 + 2 j 



(7.11) 
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Recall that on the negative ray p{S,) = Xc(Oj see (7.7). Adding (7.11) 



and (7.10) we obtain (7.6), moreover 



1 



Xc = vr(--u.(C)). 



(7.12) 

□ 



Theorem 7.1 follows from Lemma 7.2, (4.11) and (7.3). 

Completion of the proof of Theorem 1.1. 
The error term L„ satisfies 



(7.13) 



where a„ = max{0„(,^) : A < ^ < 1}. This follows from (2.4) and our explicit 
representation of the extremal polynomial (2.7). 

The necessary constants C,r],r]i,r]2, a which depend only on A and B, 
and the harmonic measure uj{x) = lo{x, [—A, — 1],C\I) were defined in the 
Introduction. 



According to (7.4), (7.12) and (7.5) we have 

X^{ll2-u:{C) + uj{D), - 
vr \-l/2 -uj{C)+u{D), - 



-K < —fl, 
-fl < —K. 



(7.14) 



Notice that w is a strictly increasing function on M\(— ^, 5) and the image 
of this set (together with the infinite point) equals [0, 1]. Therefore, for every 
n there exists a unique solution Dn of the equation 



uj{Dn) = {an + u{C)}, 

where {•} stands for the fractional part. 

Equations (4.12), (7.14) and — ^ d imply that 



(7.15) 



(7.16) 

which corresponds 



Let —Kn £ {—I, — A) be the point in C-plane (see Fig. 
to Dn in z-plane. Then g{(, — KriJ — ^ ^(C) —k). Now (7.1) implies 

hm {0n,{C) - nigiC) + g{C, -^nj) = J In ^^^^ + g{C i^)- 



(7.17) 
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The right hand side is independent of the subsequence {n/}, so the hmit as 
n — 7- oo exists in the left hand side. In the resulting formula we let (" — t- 
and obtain 



lim {9n{jy) - ng{v) + g{u, -Kn)) 



-In — +r?2. 
2 rii 



(7.18) 



The functions g{C, uniformly bounded and have bounded deriva- 

tives on (A, 1). Therefore, 

an = max{6'„(^) : A < ^ < 1} = 9n{'^) + o(l), n oo. 

Thus 

an = m- G{Dn, C) - - In + 7?2 + o(l). (7.19) 

2 rji 

To obtain the final result, this expression for an has to be substituted to 



We can simplify the expression e 



G(D„,C) 



in the resulting formula and 



avoid solving equation (7.15) in the following way. 

Let F be the conformal map of the upper half-plane onto a rectangle 
{0,p,p + where p > and the vertices of the rectangle correspond to 
(1, B, —A, —1) in this order. It is easy to see that u = Im F. So 



F{C) = iuj{C), 



and in view of (7.15) 



F{Dn) =p + iuj{Dn) =p + i{an + uj{C)]. 



(7.20) 



(7.21) 



Christoffel-Schwarz formula gives (1.2), and p = r/i, where r is defined by 
(1.4). We reflect our rectangle with respect to the imaginary axis and apply 
the map z i— )• (z7r/p)z, to obtain the new rectangle 

R = {x + iy : — 7r/p < x < 0, |y| < vr}. 

Then maps this rectangle R into a ring e^'^/^ < |tf | < 1 and we use the 
expression of the Green function of this ring [3, §55 (4)] substituting to this 
formulsj^ In If = ivr — {tt /p)u){D), Inc = {tt/p)lo{C) and using r instead of 
— 1/r. The result is simplified using Table VIII in [3] and we obtain 

^o{U{Moo)+uj{C)}-uj{C))\t) 



^9o(^({nu;(oo) + w(C)}+a;(C))|r) 



where t = ip is given by ( |1.4[ ). Combining this with (7.13) and (7.19) we 
obtain the statement of Theorem 1.1. 



^ In the English edition of 1990, this formula contains two misprints: an extra vertical 
line and missing subscript 1 in the theta-function in the denominator. 



26 



8 Example (the symmetric case) 

We consider the case / = [—A, —1]L) [1, A]. In this case 

G(z,oo) = r , =. (8.1) 



Therefore 



xdx 1 f^i^ dt 14 + 1 , , 

77=/ = - / = - In 8.2 

70 yj{x^ -l){x^ - A^) 2 A 2^-1 ' 

and 

^i = -^G"(0,oo) = ^. (8.3) 



Also, 



/-I xV(x2 - l)(a;2 - ^2) z V-42-1 



Notice that u;(oo) = 1/2, C = and uj{C) = 1/2. Therefore for 77 = 2777 + 2 

we have Dn = OO, so L2m+2 = -Z>2m+1- 

For n = 2777 + 1 we eet 



V(2777 + 1)7?. /|^e(2-+l)'^+^2L2™+l 

V 



(8.5) 



2777 + 1 + 1\'" /A+1 2A 



4A \A-lJ \l A-1 - 1 
Finally, 

hm ^2777 + 1 -^L2„+i = ./-, 8.6 

m->-oo — ly A — 1 Vvr 

as we proved in |5j . 



9 Approximation of sgn(a:) by entire functions on 

[-A-i]u[i,+oo) 

Only a minor variation of our method is needed to investigate the following 
problem: minimize 

sup{|/(a;) - sgn(x)| : x G [-A, -1] U [1, +oo)} (9.1) 
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among all entire functions f of order 1/2, type a. 



Let E{a) be the infimum (9.1). It is easy to prove the existence of an 
extremal function using normal families arguments. 

Now we describe a construction of extremal functions. We take the error 
E as an independent parameter. Let a > be the unique solution of the 
equation L{a) = E, where L{a) is defined in the beginning of Section 2. For 
h > 0, and an integer k > 2, let Hkih) = n^^(/i), that is the region in 
the upper half-plane whose boundary with respect to the upper half-plane 
consists of the segment [0, ia], and the curve 6-k as in (2.6). Let Qk,h '■ C+ — >• 



Hkih) be the conformal map normalized by @k.h{^^) = 0,Gfc,/i(oo) = oo. 

Proposition 9.1. If h = then S^QkfljO,) 'is the unique extremal function 
for 

AG[e^},{-ck),e-,},{-ck-i)]. (9.2) 

If h > then 5'(0fc^/j,a) is the unique extremal function for 

A = ek,h{-Ck-i + 0). 

The proof of this theorem is similar to the proof of Theorem 3 in [5] . We 
recall the argument for the reader's convenience. 

Proof. Let f{z) = S{@{z),a). Let xi < X2, ...—)■ +cxd be the sequence of all 
alternance points. Let u > be the type of / with respect to the order 1/2. 
Let g be an entire function of the same type a, order 1/2. Without loss of 
generality we may assume that g is real. Then there exists a sequence {yk} 
interlaced with {x^}, that is, 

xi < yi < X2 < y2 < ■ ■ ■ , 

such that f{yk) = givk)- Consider the product 

1 - z/xk 



F z) 



nr 



This product converges uniformly on compact subsets of the plane and has 
imaginary part of a constant sign in the upper half-plane and of the opposite 
sign in the lower half-plane [9l VII, Thml]. This implies that 

F(re^*) = 0(r), r ^ oo (9.3) 

uniformly with respect to t in e < t < 27r — e, for every e > 0. As f{yk) = 
giVk), we have 

fiz)-g{z) ^ Pjz) 

f'{z) {z-c)F{zy 
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where c is the critical point of / which is outside the set [—A, — 1] U [l,oo). 



If there is no such point c, then the factor [z — c) has to be omitted in (9.4). 
Then P is an entire function of order 1/2. 



Now we notice that the left hand side of (9.4) is bounded for |Im 2;| > 1. 
Indeed, g and f — g are at most of type a, order 1/2, while /' has indicator 
a"sin(t/2), < t < 2-k, so the ratio has zero type in C\M+ and thus this 
ratio is bounded by the Phragmen-Lindelof theorem. 

Combining this with (|9.3l) we conclude that P is a polynomial, and 



P{z)/{z-c) = 0{z), z^oo, (9.5) 
if the point c exists, and 

P{z) = 0(z), z^oo, (9.6) 

if the point c does not exist. 

On the other hand, P{x) = for each non-critical alternance point x. 
From our construction of / = S{Q, a) it follows that when c is present, then 
there are 3 non-critical alternance points, namely —A, —1,1, while when c 
is absent, then there are at least 2 non-critical alternance points, namely 



1, 1. Together with (9.5), (9.6) this implies that P = 0, that is, f = g. □ 



Proposition 9.2. For every E and A there exist k,h and a such that 
S{Qk,h: ^) is an extremal function for the set [—A, —1] U [1, +oo). 

Proof. For given E we can choose a such that L{a) = E. To prove the 
existence of k and h, we use a monotonicity argument as in Remarks in Sec- 
tion 3. Namely, we introduce the following order relation on the pairs {k, h): 
{k, h) -< {k' , h') if k < k' or k = k' and h > h'. With this order, the set of 
pairs {k,h) becomes isomorphic to the positive ray, and the correspondence 
(fe, h) A becomes monotone increasing. This function is continuous for 
h ^ and has a jump at each point {k, 0) (this jump is seen in the right 



hand side of (9.2)). So we can obtain any A > 1 from some pair {k,h). □ 



Theorem 9.3. For every A and a, there exists a unique extremal function 
f of type a, and f = S{Qk,h, cl) for some positive integer k, h >0 and a > 0. 

Proof. Let a{A,E) be the type (with respect to order 1/2) of the extremal 
function defined in Proposition |9.2[ Then Proposition |9.1| implies that for 
every A, the function E i— >• a{A,E) is strictly decreasing. It is easy to 
check that ct{A,1) = and a{A,0+) = +oo. Moreover, E cr{A,E) is 
continuous. So there is a unique E = E{A, a), which is the error of the best 
approximation for given A and a, and from this E and A we define k and h 



using Proposition 9.2 □ 
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To state the asymptotic result, we introduce the Martin function M{x) of 
the region C\/, where / = [—A, — 1] U [1, +c«), replacing the Green function 
which we used before. Martin's function is characterized by the properties 
that it is positive and harmonic in C\I, equals zero on / and has asymptotic 
behavior 

M{ — x) ~ ^/x, X — +CO. 

We have M{z) = Im M{z) where M is the conformal map of the upper 
half-plane onto the region 

{x + iy : x> -vra, y > 0}\[0,ir7], 

such that 

M{±1) = 0, M{-A) = --Ka, M{-x) ^ x ^ +oo. 

These relations define a and ry uniquely. 

Martin's function has a single critical point C G (—1,1) and we use the 
notation rj = M{C) and r]i = — M"(C)/2, as before. The Green function 
C) satisfies 

G{x,C) = -\n\x-C\+r]2 + 0{x-C), x^C, 

and this defines r]2- We also introduce the harmonic measure oj{z) = 
Ijj{z,[—A^—V\,C\I). Then oj{x) is continuous and strictly increasing on 
[— oo, — A), and maps this ray onto [0, 1), so the equation 

u:{D,) = {aa + u:{C)}, 

where {x} is the fractional part of x, has a unique solution for every o" > 0. 

Theorem 9.4. The error of the best uniform approximation of the function 
sgn(x) on [—A, — 1]U[1, +oo) hy entire functions of order 1/2, type a satisfies 



E{a) ^ ^{a{a))-^l\--^'^\ 



where 

aia)=w-G{D^,C)-lln'^ + m- (9-7) 

2 r]i 



The equation (9.7) is analogous to (7.19). One can simplify e'^^^'^'^^ as 
we did in Section 7 by using an expression of the Green function in terms 
of theta-functions. 
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